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Abstract 

We consider the SU(2) lattice gauge theory at finite temperature in 
dimensions, with different couplings Pt and (3s for timelike and spacelike pla- 
quettes. By using the character expansion of the Wilson action and perform- 
ing the integrals over space-like link variables, we find an effective action for 
the Polyakov loops which is exact to all orders in (5t and to the first non- 
trivial order in (3s- The critical coupling for the deconfinement transition 
is determined in the (3+1) dimensional case, by the mean field method, for 
different values of the lattice size Nt in the compactified time direction and 
of the asymmetry parameter p = ^f^tfWs- We find good agreement with 
Montecarlo simulations in the range 1 < At < 5, and good qualitative agree- 
ment in the same range with the logarithmic scaling law of QCD. Moreover 
the dependence of the results from the parameter p is in excellent agreement 
with previous theoretical predictions. 
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1 Introduction 



The aim of this paper is to obtain, by using only analytical methods, reliable esti- 
mates of the deconfinement temperature in the SU(2) pure gauge theory (namely 
without quarks) in (3+1) dimensions. The natural framework to pose this question 
is that of the finite temperature Lattice Gauge Theories (LGT). In this framework, 
during these last years, the best estimates of the deconfinement temperature have 
been obtained by means of Montecarlo simulations, which are certainly the most 
powerful tool to extract quantitative results from LGT. However we think that it 
is important in itself to have some independent analytical estimate of the location 
of the critical point, besides the outputs of the computer simulations, to reach a 
deeper theoretical understanding of the deconfinement transition. The attempts 
to obtain analytically the critical temperature have a rather long history, starting 
more than ten years ago ||, |^, H, |^ . However the strategy has always been essen- 
tially the same: first, construct an effective action in terms of the Polyakov loops 
(which, as we shall see below, are the relevant dynamical variables in the physics 
of deconfinement for pure gauge theories). Second, use a mean field approximation 
to extract the critical coupling. A common feature of all these attempts was that 
the effective actions were always constructed neglecting the spacelike part of the 
action. As a consequence it was impossible to reach a consistent continuum limit 
for the critical temperature. 

The aim of this paper is to show that it is possible to overcome this problem. 
We shall construct in the SU{2) case an improved effective action which takes into 
account also the spacelike part of the original Wilson action and is exact to all 
orders in the timelike coupling. This is a rather non trivial result and we shall 
devote most of this paper to describe how it can be obtained. Moreover, as we 
shall see, our approach is a constructive one and can be extended in principle to all 
orders in the space-like couplings. 

We decided in this paper to concentrate only on the gauge group SU{2) for 
simplicity reasons, but most of our results can be extended to SU {N) models with 
N > 2. Indeed, this paper can be considered as the natural continuation of |^ 
where these same techniques were applied to the N ^ oo limit of LGT. Here we 
try to eliminate the large N approximation by looking directly at the N = 2 case. 

This paper is organized as follows: after a short introduction to finite temper- 
ature lattice gauge theory (sect. 2), we shall devote sect. 3 to the construction of 
the effective action. In particular, sections 3.3.1 and j.3.2 contain the computation 
of the first non-trivial contributions from the space-like part of the action; these 
sections are rather technical and the reader interested mainly in the results may 
wish to skip them, as the results are anyhow summarized in section p.4| . In sect. 4 
we shall extract the critical deconfinement temperature with mean field techniques, 
we shall discuss our results in comparison with existing Montecarlo estimates and 
check their consistency in the case of asymmetric lattices with known theoretical 
results. Finally sect. 5 will be devoted to some concluding remarks. We shall try 
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to keep our formalism as general as possible, so we shall derive in sect. 2 and 3 
the effective action for the Polyakov loop in a (rf + l)-dimensional LGT with an 
arbitrary and we shall fix = 3 only in sect. 4. 



2 Finite Temperature LGT 
2.1 General Setting 

Let us consider a pure gauge theory with gauge group SU{2), defined on a + 1 
dimensional cubic lattice. In order to describe a finite temperature LGT, we have to 
impose periodic boundary conditions in one direction (which we shall call from now 
on "time-like" direction), while the boundary conditions in the other d direction 
(which we shall call "space-like") can be chosen freely. We take a lattice of Nt (Ng) 
spacings in the time (space) direction, and we work with the pure gauge theory, 
containing only gauge fields described by the link variables Un-i G SU{2), where 
n = (x, t) denotes the space-time position of the link and i its direction. It is useful 
to choose different bare couplings in the time and space directions. Let us call them 
Pt and Ps respectively. The Wilson action is then 



Sw = T.l{(^tYl Trf(t/„;o.)+/3sE Trf(f/„,,, 
^ I ^ i<j 



where Trf denotes the trace in the fundamental representation and Un-oi {Un;ij) are 
the time-like (space-like) plaquette variables, defined as usual by 

Un]ij = Un-^iUn+rjU.^_^_j.jU^.j . (2) 

In the following we shall call Sg (S-t) the space-like (time-like) part of Sw 
Let us introduce an asymmetry parameter p defined by the relation: Pt/Ps = P^- 
As p varies we have different, but equivalent, lattice regularization of the same 
model. This equivalence is summarized by the following equations, which can be 
obtained by taking the classical continuum limit of (^ and which relate Ps and Pt 
to the (bare) gauge coupling g and to the temperature T: 
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2 a^-yPsPt , T = --^l^ . (3) 

Here a is the space-like lattice spacing and is the time-like spacing, hence p 
is the ratio between the two. From this last observation it is clear that equivalent 
regularizations with different values of p require different values of Nt. Hence, to 
maintain the equivalence, Nt must be a function of p: Nt{p). 

Among all these equivalent regularizations a particular role is played by the 
symmetric one, which is defined by: 

P = -.a"-' (4) 
9^ 
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(from now on we shall distinguish the symmetric regularization from the asymmetric 
ones by eliminating the subscripts t and s in Comparing eqs.(||,|) we see that 
all the regularizations are equivalent if the following relations hold: 

P = pf3s = ^, (5) 
P 

N,{p)=pN,{p=l). (6) 

Notice however that these equivalence relations are constructed in the naive or 
"classical" continuum limit. At the quantum level, in the (3+1) dimensional case 
these relations change slightly. These modifications have been studied by F.Karsch 
in 1^ and we shall discuss them in sect. 4 . Let us now only anticipate that eq.(|D 
becomes: 

A = p(i3 + 4c,(p)) (7) 
A = ^ii^^, (8) 

P 

where the two functions Co-(p) and Cr{p) can be found in |^ []. For our purposes we 
only need to know the first terms of their 1/p expansion in the p — >■ cxd limit. Let 
us define: 

4c.,. = <. + ^ + ---. (9) 

The a's can be calculated from the expression of c^^rip) given in 0; they are: 
a° = -0.27192; al = 1/2; = 0.39832; al = 0. 

These relations will play a major role in the following, since the invariance of our 
results as p changes is a crucial consistency check of all our approach. Actually it is 
a non trivial test, since these relations are the result of a weak coupling calculation, 
and only become manifest in the continuum limit of the model. Being able to 
reproduce them within the framework of a strong coupling calculation would be a 
remarkable and a priori unexpected result. This is actually the discussed 
in detail in sect. 4. 

A second reason for which it is important to have under control this p symmetry 
is that, at the end, we would like to compare our prediction with the Montecarlo 
simulations, which are all made on symmetric lattices. However, at the same time, 
it is only in the limit of highly asymmetric lattices (p oo) that (as we shall see 
in more detail below) we can define in an unambiguous way our expansion of the 
spacelike part of the action, so it is somehow mandatory for us to be able to match 
these two different limits. 

In a finite temperature discretization it is possible to define gauge invariant 
observables which are topologically non-trivial, as a consequence of the periodic 



^We have chosen, for sake of clarity, to keep the notation of for and Cr- This explains 
the factor 4 in the above equations 
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boundary conditions in the time directions. The simplest choice is the Polyakov 
loop defined in terms of link variables as: 

Nt 

P^ = TTl[{V^,t) . (10) 

t=i 

where V^^t = Ug^ffi are the vertical link matrices. In the following we shall often use 
the untraced quantity P^, defined as: 

Nt 

Ps=l[iV,,t) . (11) 

t=l 

which will be referred to as "Polyakov line" . 

As it is well known, the finite temperature theory has a new global symmetry 
(unrelated to the gauge symmetry), with symmetry group the center C of the gauge 
group (in our case Z2). The Polyakov loop is a natural order parameter for this 
symmetry. 

In > 1, finite temperature gauge theories admit a deconfinement transition at 
T = Tc, separating the high temperature, deconfined, phase (T > Tc) from the low 
temperature, confining domain (T < T^). In the following we shall be interested in 
the phase diagram of the model as a function of T, and we shall make some attempt 
to locate the critical point T^. The high temperature regime is characterized by the 
breaking of the global symmetry with respect to the center of the group. In this 
phase the Polyakov loop has a non-zero expectation value, and it is an element of 
the center of the gauge group (see for instance 0). 

2.2 Svetitsky-Yaffe conjecture 

The Svetitsky-Yaffe conjecture is based on the idea that, if one were able to 
integrate out all the gauge degrees of freedom of the original {d + l)-dimensional 
model except those related to the Polyakov loops, then the resulting effective the- 
ory for the Polyakov loops would be a dimensional spin system with symmetry 
group C. The deconfinement transition of the original model would become the 
order-disorder transition of the effective spin system. This effective theory would 
obviously have very complicate interactions, but Svetitsky and Yaffe were able to 
argue that all these interactions should be short ranged. As a consequence, if the 
transition point of the effective spin system is of second order, near this critical 
point, where the correlation length becomes infinite, the precise form of the short 
ranged interactions should not be important, and the universality class of the de- 
confinement transition should coincide with that of the simple spin model with only 
nearest neighbour interactions and the same global symmetry group. In particular 
the deconfinement transition of the d + 1 dimensional SU(2) LGT in which we are 
interested should belong to the same universality class of the magnetization tran- 
sition of the d- dimensional spin Ising model. Unfortunately this argument cannot 
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help to fix the critical temperature, which is not an universal result, but depends 
on the precise form of the action that we study, and hence of the short ranged 
interactions that we neglected above. In the next section we shall construct these 
correction terms explicitly. 



2.3 Character expansion 

An important role in the following analysis will be played by the character expan- 
sion, which in the SU(2) case is very easy to handle. Let us briefly summarize few 
results. The character of the group element U in the j*'* representation is: 

where Trj denotes the trace in the j*'* representation and 6 is defined according to 
the following parametrization of U in the fundamental representation: 

f/ = cos(^)l + ia-nsin(^). (13) 

Here f? is a tridimensional unit vector and cTj are the three Pauli matrices. Notice, 
as a side remark, that with this parametrization the Haar measure has the following 
form: 

Z^f/ = sin2(^)__ (14) 

and the Polyakov loop becomes = 2 cos{6£) 

The following orthogonality relations between characters hold: 

j DUxriU) XsiU) = 5r,s (15) 



J2drXriUV-') = 5iU,V) (16) 

r 

where dr denotes the dimensions of the r*'* representation: = 2r + 1. In the 
following we shall use two important properties of the characters: 

/ DUxriU) Xs{U~'V) = Sr,s^^ (17) 

J DUXr{UV,U-W2) = jXr{V,)Xr{V2). (18) 

The character expansion of the Wilson action has a particularly simple form: 

ef^^(^) = X:2(2, + l)^^X.(t/), J = 0,^,1--- (19) 

where /„(/?) is the n^^ modified Bessel function. It is customary to collect in front 
of expression ([T9|) a factor of "^^^/f^, so that the expansion starts with 1. 
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3 Construction of the Effective Action 



In this section our goal is to construct an effective action for tlie finite temperature 
LGT in terms of tfie Polyakov loops only. To do so one should be able to integrate 
exactly on the spacelike variables so that the only remaining degrees of freedom 
at the end are the Polyakov loops. Notice that in this way the resulting effective 
action would live in d dimensions (one dimension less than the starting model). 
This is exactly along the line of the original Svetitsky-Yaffe program. As already 
remarked in the introduction the early attempts to determine analytically the criti- 
cal temperature were all based on the assumption that the deconfinement transition 
is dominated by the timelike plaquettes, and the contribution of the space-like pla- 
quettes was consequently neglected. Although this approximation correctly predicts 
the existence of the deconfinement transition, a quantitative estimate of the crit- 
ical temperature for large enough values of iVj, namely near the continuum limit, 
requires the contribution of the space- like plaquettes to be taken into account. Ac- 
cordingly, we shall treat the timelike part of the Wilson action St as a Born term 
and treat the spacelike part Ss as a perturbation; namely, we shall make a strong 
coupling expansion in f3s-, while the time-like part of action will be treated exactly. 
This means that order by order in j3s the dependence of the effective action from 
(5t will be exact, the only expansion parameter being thus (5s. Of course, the zeroth 
order in fig will contain the timelike plaquettes only. It is not at all obvious that the 
integration over the spacelike links could be done to all orders in but it turns 
out to be the case in the framework of the characters expansion (sec below) order 
by order in /3s- Rather than a straightforward expansion in powers of (3s we shall 
use for each space-like plaquette a character expansion. Each representation j in 
the expansion gives a contribution proportional to a ratio of Besscl functions which 
is of order so that the character expansion and the expansion in powers of (ig 
coincide up to higher order terms arising from the power series expansion of the 
Bessel functions. As it will be discussed later in sect. 4 these higher order terms 
vanish anyway in the limit of highly asymmetric lattice. We shall consider in this 
paper only the zeroth order and the first non trivial order in (3s, namely Pg- Terms 
of order (3s come either from one space-like plaquette in the adjoint representation 
or from a couple of plaquettes in the fundamental representation. However it should 
be noted that there is no obstruction in principle to go to higher orders. 

For any given order in (3s the result is given by an infinite sum of characters. 
Remarkably enough in the Nt = 1 case this series can be summed exactly and the 
result can be written in a closed form. This is essentially due to the fact the if 
Nt = 1 the same effective action can be obtained in a completely different way, 
using techniques typical of matrix models (see below), thus allowing a non trivial 
check of all our strong coupling results. 
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3.1 Expansion in f3s of the effective action 

The effective action Ses for tlie Polyakov lines = Hf^i is obtained by inte- 
grating over all the spacelike degrees of freedom in the action (|ID. As explained 
previously, our approach is to consider the contributions from the spacelike plaque- 
ttes up to a certain order in Ps only. So, for our purposes, it will be convenient to 
expand separately the spacelike and the timelike part of the action (|1]): 



exp (5, 



cffJ 



/ n DUs^t;i exp Sw 

n Du^',t'■,^' n 



x't':i' 



x".t";i" 



1 + ^ dj^YJ^Y^Xj{Ux",t";Oi"] 




Xl{Ux,t; 



I]) 



(20) 



Specifically, we work out here the effective action up to 0{l3l). This means that in 
eq.(P(]D we must look only at the terms containing at most a single space-like pla- 
quette in the adjoint representation, Xi{Us,t;ij), or two space-like plaquettes in the 
fundamental, Xi{U£,ti;ij)x^{Uy,t2;ki)- -Due to the orthogonality relations for char- 
acters, it easy to convince oneself that a pair of plaquettes in the fundamental 
representation do actually contribute to the integral only if they appear in the 
same spatial position (at two different times ti and ^2); for the same reason a sin- 
gle fundamental plaquette cannot contribute. We are thus lead to the following 
expression: 

exp(s'eff) = / n ^f^w n fi+f:rf /';7iy ^^(^^"'.^";o^"> 

^ TV, . / ^ N / . / ^ N \ 2 



x,i<j 



t=i -'1 



I'M 



-XiiUx-., 



t;ij. 



+ E4 

tl<t2 



J 



MPs 



Xl{U£,t^;ij)Xl{U£,t2;ij) 



(21) 

In the next sections we will consider separately the three contributions appearing 
in (PTf)- The first one that we will consider corresponds to the "1" in the second 
factor above, and gives the 0(/?°) result. 



3.2 Zeroth order approximation 

In the zeroth order approximation we have to consider only the timelike part of the 
Wilson action: 



exp(S'o) 



n 

x,i:t 



Da 



x,i;y 



\ 



(22) 
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In this case we can easily integrate all the spacelike links. The reason is that 
each spacelike link only belongs to two timelike plaquettes; hence, by making a 
character expansion, it can be exactly integrated out. Let us do this integration in 
two steps, for future commodity. First let us integrate all the spacelike links except 
the lowermost ones (which, due to the periodic boundary conditions coincide with 
the uppermost ones). We obtain, using eg. (|T7|) : 



exp(5o) = n 1 + E 



hj+i{Pt 



\ 



(23) 



where is the open Polyakov line (whose trace is the Polyakov loop) in the site x 
and Uss-i are the remaining lowermost spacelike links. Integrating also on U^-i using 
eq.(|T8D we end up with 



exp(^o) = n 1 + E 



hj+i {Pt 



\ 



(24) 



Let us define, for future convenience, the linkQ element of exp(5'o) as follows: 



j=Q L 



Nt 



(25) 



It is now evident that this basic element, which will be denoted also as = 
C^{Oxi 6x+i), depends only on 6^, 6g+i, which are the invariant angles for the Polyakov 
lines Pg, Pg+i in the sites joined by the link. Indeed from now on we will always 
assume to have gauge-rotated the Polyakov lines to be diagonal: 







(26) 



Notice that the zero^'^-order action (^4]) is simply given by 



(27) 



3.3 First order approximation 

The 0{(3^) effective action (|27|) contains just nearest-neighbour interactions be- 
tween the Polyakov loops. As we shall see in the following, the net outcome to the 
effective action from the 0{j3g) terms in ( pT]) is the addition of interaction terms 
involving more than two Polyakov loops. Specifically, the term of order /3g contains 

^The links we are referring to are those of the d-dimensional spatial lattice, corresponding to 
a space-like slice in the original d + 1-dimensional lattice 
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single plaquette in the 
adjoint representation 



pair of plaquettes in the 
fundamental representation 




interactions among the invariant angles of all the Polyakov lines around a spatial 
plaquette. 

Let us now compute the contributions from the adjoint space-like plaquettes and 
from the pairs of fundamental ones in eq. (pT]). As already said in the introduction, 
the next two subsections, containing these computations, are quite technical; the 
results are summarized in section p.4| . 

3.3.1 The adjoint representation term 

To calculate the contribution from the adjoint representation, we have to select 
from (|2l|) the term: 

h[Ps) -I ^/ j/y x",t";i" \^ j=i ^^yPi) j 

Nt 

X E T.Xi{U,,;,). (28) 

x,i<j t=l 

Let us fix a spatial position x,i < j in the above sum over the space-like plaquettes, 
and study the corresponding integral. The integration over all the link matrices not 
pertaining to the chosen spatial position can be performed exactly as in the 0{f3g) 
case, giving as a result the product of all the link factors (7°, ■, except those in the 
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chosen spatial position^ (i.e. except C° ., C°_^jj, C^_^_i_^_j_i,C^_^_j _j). To treat the 
remaining non-trivial integrations, first we can note that all the spacelike plaquettes 
in the same spatial position give evidently the same contribution, regardless of the 
time t; therefore the sum over the time positions in (^81) results simply in a Nt factor. 
Secondly, it is convenient to use the following relation for the SU{2) characters: 



Xi = (Xi)' 



(29) 



The "—1" simply reproduces the zeroth order term, and gives a renormalization of 
order (3^ to such contribution. The integral over the link variables along the pla- 
quette can now be decoupled into products of integrals over single link matrices, by 
writing explicitly [xi{Ux,t;ij)]'^ in term of traces of link variables in the fundamental 
representation. Thus eq.(^) can be rewritten in terms of the following integrals 
over the unitary spacelike link matrix U: 



B. 



(p^,p^+,) = Jdu ii + Y.d, 



Nt 



XAUPx+^U^Pl:) U^pUls (30) 



where a, . . . = 1,2 are the indices of the U matrix in the fundamental representation. 
Let us assume that the Polyakov lines in eq.(|30D are already set in the diagonal 
form of eq.(^). Then the measure and the argument of Xj at the r.h.s. of eq.( 
are invariant under the transformations 



and 



Uaf3 — ^ ^aa Uaf3 , 
Ual3 — > Ual3 , 



(31) 
(32) 



'-fS ~~' )ll ^yS 

where is a diagonal SU (2) matrix. By using this invariance one can easily conclude 
that Bafjys = unless (3 = j and 6 = a. As a consequence, the integral ( PU[ ) depends 
on the invariant angles of the Polyakov line only, and can be written as follows: 



Bap-ysi^Px^ Px+'l) = Ba/Sysi^xj ^x+i] 



/3-yOSa '-^a/Sy^xj ^x+i 



(33) 



(no summation over repeated indices). Moreover, it is not difficult to show that 
Cai3 is a real symmetric matrix. By using these facts, we can write the contribution 
(P^) to the effective action in terms of the invariant angles of the Polyakov lines. 
The contribution to eq. (|28|) at fixed x, i can be expressed as: 



3iVy 



hi/3s 



(34) 



^Later on we will denote this product as Ilimk^pi C'o(link). 
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where the term (-1) in ( |3^ corresponds to the term (-1) in (|29| ) and the matrices 
C{6^, Ox+i) are a normahzed version of C: 

^a,l3[t'x,fj!+i) — — Q {OD) 

Indeed in writing eq.(p^) we have muhiphed and divided by C2^- C°+ij- C*°+j+j,_i- 
C^j^j _j in order to collect the factor nx,jC*9,j = exp^o- 

The last step is the explicit evaluation of the matrix elements Cap. We set 



j=0 



Nt 



C% , (36) 



so that C^^p is defined as the contribution of the f"^ representation in eg. (pO|) : 

I DUUafiU\,xA^Ps^,V^Pt) = SasS.pC^l (37) 

It follows from (|37D that the non vanishing integrals at the l.h.s. depend only on 
|[/„/3p, and hence that C{{^ = C^^^ and CjJ = C^{\ 

To compute the matrix elements C^^^ , which is not a completely trivial task, we 
use the following strategy^. We note that the matrix C^^'^ can be expressed in terms 
of the integral 

K^'KOx. 0,+.) = I DUxj{UP2U^Pl) = jxAP2)xAPl) (38) 

through a system of two linear Schwinger-Dyson-like equations. Indeed, considering 
the integral / DUUap Up^ Xj{U Ps+iU'^ pI) we easily find that 

Cj? + Cjj = ir(^'\ (39) 
To construct a second independent equation, let us consider the integral 

j DUxL^iUP,+,U^Pl)xjiUP,+,U^Pl). (40) 

On one hand we can write the character xi explicitly as a trace and express the 

integral in terms of the C^p by using eq. (pT]) . On the other hand, the integral (|40|) 
can be written in terms of K^^^ functions by using the basic SU(2) Clebsch-Gordan 
relation: X^Xj = Xj+i + Xj-i- The resulting equation is: 

2 cos(^,^+, - e,)C^^ + 2 cos(^^,^+, + e,)C^^ = K^^-"^) + K^^H'^ (41) 



^There is another possible way to compute C^^^, based on the expression of the SU(2) characters 
as Tchebicheff polynomials of second kind: Xj(0) = U2j (cos(6')), which was utilized in However 
this alternative technique cannot easily be applied to the case of a pair of fundamental plaquettes. 
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Eq.s ([39| ) and (^) form a set of two linear equations in the two unknowns Cii and 
whose solution is: 



(i) 



'-'12 



4 sin Ogj^i sin 6'^ 

/sfO-l) - 2cos(^£+i - + 

4 sin sin 6'^ 



(42) 



By inserting these results in eg . (^61) we finally obtain the coefficients. We 
choose to write the matrix C in the form 



^ f Cq + Ci Cq — Ci 
2\Cq — Ci Cq + Ci 



(43) 



which is consistent with the symmetries of C and allows an easy evaluation of the 
trace contained in eq. (|3^). After some algebraic rearrangements it follows from 
eq.(|D that 



1 



2 sin sin d,i^ 



E 



(2j + l)/i(A)^' 



Nt 



+2cos[(2j + l)fe+i] cos[(2j + l) 
+2 cos 6*^+4 cos 6*^ - 



+ 



T \ 



h 



(44) 



while 6*0(6*^+4, 6'^) coincides with the link element as defined in eq.i 

With the parametrization of the matrix C given in eg. (^3]) the trace in eq.i 
simply reads: 

= Co{9g, 9£+i)Co{9g-^i, 9g^i-^-j)Co{9j;^i^j, 9g-\^j)Co{9rjt^j, 9^) 

+ Cl(6'i', ^:f+j)Ci (0^4.4, ^:r+i+i)Cl(^x+i+j) ^x+j )C'l (^x+j ) ) • (45) 



By inserting this expression into eg. (|3^) , the products of Cq's cancel, and we are 
left with the explicit expression for the contribution of an adjoint plaguette: 



hiPs 



where 



11^^,4 Ci{9g, ^:f+i)C*l(^x+i) ^x+i+j)C*l(^x+i+j) ^x+j)C*l (^x+jj 9g) (46) 

Ci{9^, 9g^i) 



Ci{9g, 9^j^i 
and Ci{9^, 9^+i) is given by eg.(|4|). 



(47) 
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3.3.2 Pair of fundamental representations 

Let us go back to eg . (|2l|) , and consider the last type of contributions, namely the 
ones coming from two plaquettes in the fundamental representation. As in the 
previous case of the adjoint plaquettes, we consider the contribution of a single pair 
of plaquettes; that is, we fix the spatial position x,i < j of the plaquettes. We 
can moreover fix the time position of one of these two plaquettes, say at ti = 0. 
Then the sum over ti just gives a factor of Nt. The second spatial plaquette will 
be located at t2 = M, (M = 1, . . . , Nt — 1). Notice that the contribution from the 
fundamental plaquettes is not present when Nt = 1. 

To perform the computation of this contribution, it is convenient to choose the 
gauge so that the Polyakov lines are concentrated for instance in the uppermost 
vertical links; this choice is always possible. 

The integration of all the space-like links not involved in the two space-like 
plaquettes can be performed in the usual way. With notations similar to those of 
eq.(pOD [see also Fig. 1] we are left with the expression: 



n cdiuk) X 4(M|))^iv, 

link^pl -^IKHs) 



X 



M=l x^yepl 



1+ E 



'2/+1 



(A) 



Nt-M 



"MAT 



M 



XmiUx,yU,ji^g, 



(4J 



By writing explicitly the products in the traces, one can express eq. 
integrals of the form 



in terms of 



(49) 



where the Kronecker deltas at the r.h.s. originate, as in the case of the adjoint 
representation, by the symmetry given in eq.s (^[321) . We have assumed in eq. (|49D 
to have diagonalized the Polyakov lines as in eq. PB|) . The original expression 
can then be written as 



n Co(link)x4iV,(^) 



link^pl 



X 



M 



(50) 
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where the 2x2 matrices C^'^^ are given by 



-^2m+l(A) 


M 


-^2n+l(A) 


Nt 


[ im \ 




[ im \ 





c 



{m,n) 
a/3 • 



(51) 



Unhke the C^-jJ defined in eq.(^), C^™'"'* is not a real symmetric matrix. It rather 
satisfies the following properties: 



'-11 



'-'22 



' p{m,n) 
'-'12 



'-'21 



Moreover, one can easily show from its definition (|49|) that 



' p(m,n) 



•'a/3 



(52) 



(53) 



andCi7V.,^.+.)=4r^( 

In order to compute the matrix elements of C^"^'"''' we follow a strategy analogous 
to that utilized in the adjoint plaquette case. We show how the Cj^'^\ that have 
four independent real components, can be expressed in term of the integral 



6. 



5r, 



DUDU x^{UU^)xn{UPxMU^Pl) = -^Xn{e,+i)Xn{e,) ^ (54) 



by means of a set of four linear Schwinger-Dyson-like equations. Consider first the 
integral 

" DU DUx-L{UU^)Xm{UU^)Xn{UPx+^U^Pl). (55) 



On one hand we can write explicitly as a trace the xi factor and use the definition 
(^); on the other we can use the Clebsch-Gordan relation to rewrite the integral 
(p5|) as a combination of K^'^'> functions via eq.(0). We thus find: 



^(m,n) ^(m,n) ^(m,n) ^(m,n) _ / r n 1 ^(n) 

'-'11 +'-'12 +'-'22 +'-'21 — lOm+in + "^m-inJ 



(56) 



By considering an integral analogous to (0), but containing the trace xi (UP^-^-iWP^) 
instead of xiiUU"^), we obtain that 



(57) 



To determine all the components of C^™'"^ we need two more relations, that can be 



15 



obtained by taking derivatives of eq.(|5^) with respect to the invariant anglesQ: 

i{2n + l)[e'('^^'^^^^cir'> - e-'(^^^+^^^-)c{-''^)] = ^{5^,^^ - S^,n-0d^K^-^ (58) 
and 



i(2n+l)[e^(^--^^^)c;r'"^-e-'(^--^^^C(r^] = ^{5^,r.^ - S^,n-0d.K(-\ (59) 



where d± stands for 0/ ± We can now obtain the expression of C^'"^'' by 



(m,7i) 



de. 



a/3 



solving the system formed by the four equations (|56|j57|j58|j59| ) . This is more easily 
done in terms of the matrix C^™'"-), defined by = e'^^^+^'^^^C^™'"^ and cj^'"'' = 

e-i{e2+es+.)cj^'"\ i.e. by 



p(m,n) _ f p \ f p \ p{m,n) 



■-0/3 



(60) 



The matrix C^'^'"'' enjoys the same symmetries (^) as Cj^'"'\ while eq. (|53|) is re- 



placed by 



(61) 



The symmetries (|5^) can be implemented by writing the 2x2 matrix as 



1 I /7(™->") -KjC*","-) /^('Tij") I -KjCm,") 

~ ' -i{m,n) p(m,n) ■ Y^{m,n) p{m,,ri) . p(m,n) . ■ yy{m,n) I ' / 



2 y Cq — Ci — Cq + Ci + iB 



The equations (|56|j57|j58|j59|) imply that the only non-vanishing components of C^™'"-* 
are those with n = m ± |. The solution of the system can be expressed as follows: 



2d„ 



2(2m + 2) 



(m,rra— i) 



2(2m; 



(m,m+^) COS [(2m + l)^^+j] COS [(2m + 1)^^] + cos^^+j cosO^K^"^^ — 

(2m + 2)2 sin sin 9^ 
.1) cos [(2m + l)^s+i] cos [(2m + 1)^^] - cos^^^+i cosOgK^""^ + K^"'-^^ 



(m,m- 



(2m) 2 sin sin 6g 



(63) 



^SU(2) characters satisfy [xj4.i(^^) — Xj-i(^)] = (2j + l)xj(^), as can be seen by 

2 

expressing them as Tchebicheff polynomials: Xj(^) = U2jicos{9)). One has then for instance: 
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In analogy to eq.(^TD we can introduce a matrix defined by the relation 



-^2m+l(A) 


M 


-^2n+l(A) 


Nt 


[ hm \ 




[ A (A) J 





We can now insert the matrix C*^*^^ instead of C''^^^ in the expression (|50D , as 
the extra phases appearing in C^^^") cancel in the trace. 



The matrix C^,^ can be written, in analogy with eq. 



as 



q{M) 



.(M) ^(Af) _ , r(M) 



(M) , j^(A/) 



The explicit form of Cr' and 5^"^ can now be obtained from eq.s 
(|63). The result is: 



(65) 
and 



(M) 



1 OO r r T ^ 

i 1 ^ -19 



E 



+ 



/l 

1 oo 

- y 



'2m+l 

Nt-M 



'2m+2 



Nt-M 



+ 4 







7Vt-M\ 




-^2m 




1 






2 







-^2m+l 


M / 


-^2m+2 




-^2m 


Nt-M\ 


i h \ 








i h \ 





'x+ii ^x ) 



(M) 



1 



2 sin Ojt^i sin 9^ 

CO 

+ E 



"/2- 


M 






+ 






Jl. 



Nt-M 



-^2m+l 


Nt-M / 


-^2m+2 


M 


l2m 


M\ 


h \ 




[ h \ 




V h \ 





+ 
+ 



'2m+l 



h 



'2m+l 



h 



M 



M 



l2m-\ 


-2 


Nt 


-M 


l2m 


h 






+ 


V h \ 


-^2m4 


^2 


Nt 


-M 




h J 






i h \ 



2 cos 9^^i cos 9^ 

cos [(2m + l)6'^+j] cos [(2m + 1)9^ 
cos9,+,cos9, K^^\9,+,,9,] 



Nt-M-^ 



(66) 



In the above formula, the argument of all the Bessel functions is the timelike- 
coupling Pf 

The explicit expression (|66|) of the matrix C^^'^^ is quite complicated; but there 
are some non-trivial consistency checks that we can perform on it. First, it is not 
difficult to see that if we set M = in the formulae (|66|), namely if we let the two 
fundamental plaquettes coincide at the same time position, we correctly reproduce 
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the results obtained in section |3.3.1| for the xiXi term in the adjoint plaquette: 



%i (^x+i) ^ 0. 



(67) 



Second, we can rearrange, through some algebraic manipulations, the expressions 
in such a way that the symmetry under the exchange of M with Nt — M : 



.Ox+i)]* - iPg+i)f3l3{Px)aaC^' ^^\6g,6£^ 



(6J 



which comes from the analogous property ( ]61| ) of Cj^'"'\ becomes manifest. It is 
possible indeed to write C*-^^-* as a matrix in the following form: 



gi{2m+l)6»_ _g-i(2m+l)e+ \ 


l2m+l 


M 




Nt-M 


_gi(2m+l)e+ g-i{2m+l)e_ ) 


[ h \ 




i h \ 





4 sin 9^ sin 

^ f,i{2m+l) 

+ 



E 



-i(2m+l)0_ 



,i(2m+l)04 



_g-i(2m+l)6»+ gi(2m+l)6»_ 



'2m+l 



M 



'2m+2 



Nt-M 



1 




-^2m+l 


Nt-M / 


-^2m+2 


M 


l2m 


M \ 


-1 




V h \ 




i /l J 




V h \ 





AO. 



— e" 



-i04 



+ 



1 



-1 



-i6»_ 



-^2m+l 


M / 




Nt-M 




Nt-M 


V h \ 




V h \ 




[ h \ 





-1 1 



h 



M 



+ 



COS 

- COS 6 



+ 



- cos 9^ 
cos 9_ 



-I2 


Nt-M \ 


Jl. 


i 



(69) 



Finally, having determined the matrix C^^'^\9^+i, 9g) we can insert its expression 
into eq.(|50|), as we discussed above, to get the contribution of a pair of fundamental 
plaquettes to the effective action. The result of the trace is cumbersome and not 
very illuminating, so we shall not write it here; it will however be used in the mean 
field analysis of the following sections. 



3.4 The effective action up to 0{j3l) 

Let us summarize here our results by reporting the form of the effective action for 
the Polyakov loops determined in the previous sections. To this action we will in 
the later sections apply standard and improved mean field techniques in order to 
extract the value of the critical coupling. We have: 

exp Scflf = exp(S'o + ^i) (70) 

where 

exp So = Y^ Co{9g, 9s+i) (71) 
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and 

Si 




'XI "X+l 



(72) 



The quantities in the above equations are defined as follows: CQ{9:g,9g^i) = C9j is 
given by eq. (^5]) , Ci{9^, 9g+i) by eq.(^lD and the matrix C^^\9^, 9^+i) by equations 
(|65|) and ( ]66| ) or, alternatively, by eq. (|69|) . 

5*0 is the 0(/3°) effective action. It contains a sum over the links x,i in the d- 
dimensional spatial lattice and each term of the sum represents a nearest neighbour 
interaction between Polyakov loops. 5*1 describes the effect at 0(/5f ) of the space-like 
plaquettes, and it is given by a sum over the plaquettes x,i < j in the d-dimensional 
lattice. Each term represents the interaction among the four Polyakov loops at the 
vertices of each plaquette. Notice that in eqs. (py-|72|) the contribution of the 
space-like plaquettes has been exponentiated, which is correct at the order 0{P]) . 



3.5 The Nt = 1 case and the Kazakov-Migdal model 

The interesting feature of the A^t = 1 case is that the model that we are studying 
becomes a particular case of the Kazakov-Migdal model [0. This connection was 
already noticed in and was the origin of our previous analysis in the N ^ oo 
limit 0. All the integrals that we have described in the previous sections can 
be directly evaluated in the Nt = 1 case as particular instances of a nontrivial 
generalization of the so called Itzykson-Zuber integral |]13[, evaluated in []T^. This 



alternative derivation in the A^^t = 1 case provides another check of our computations 
(at least for the contribution of the adjoint plaquettes, as this is the only 0(/?g) 
contribution when A^^ = 1). 



The basic integral used in the computation of Sq [see eq.(|2^)] coincides, when 
Nt = 1, with the link integral 



aft cos(92-es+i) _ pft cos(0s+0s+j) 

(73) 



2(3t?,m.{9s) sin(6'^+i) 

which was non-perturbatively computed in We can compare this expression 
with our general result (P^. We find that, if A't = 1, the character expansion 
contained in eq.(p4D can be summed exactly. In fact by inserting the explicit form 
(|l^) of the characters into eq. (p4|) and then using the relation 

2 sin[(2r + l)^s] sin[(2r + l)fe+,] = cos[(2r + 1)(^,^ - 0,^+,)] - cos[(2r + 1)(^,^ + ^s+i)] 

(74) 
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and the well known expansion 



cos i 



Ioif3) + 2j2lki(3)cosike) , (75) 

k=l 



it is easy to obtain: 



g/3t cos(eg-e3+j) _ ^i3tcos(9g+es+i) 

exp(S'o) = n — AT ff^\ ■ fo \ ■ (a — ^ — • ^'^^ 
-l-^ 4/1 ^m.{Bi) sin(6'^+i) 

This expression coincides with eq. ( |75D , except for the irrelevant overall factor "^^^ff^ 
[see the remark after eq. (|T9|) ]. 

In the Nt = ^ case, the first non-trivial contributions from the space-like pla- 
quettes, i.e. those coming from an adjoint plaquette, can be extracted from the 
definition of the correlators given in 



JdU,.,exp {f Tr; {P,U,;PV^IM;)} (U^.) {uj 
J dU,, exp {f Tr^ (P^f/.^.Pt+.f/lJ } 



p,a 

— (77) 



These correlators correspond to the B^i^p„ of eq. (|30|) , divided by C^ f, they must be 
therefore diagonal, namely: 

{iU^.i),,u = ^lKC,,,{x- 1) (78) 

where the Cp^^p{x\ i) are equivalent, apart from the different normalization, to our 



Chi matrix elements. These correlators were calculated in [14 



2l3t sin(e^) sin{e^+i) - ( 1 - e'^^* ''^'^^-^ sin{e2+,) 



(^l _ e"2A«i°(es)sin(ej+,)^ {2(3tsm{9g) sin(e^.+i)) 



A A .^ 1 - e-2/^'«'"('^^-)^i"(^^-+^) (1 + sin(fe) sin(^,.+,)) 

Ci2(x\i) = C2i[x;t) = —p r (79) 

' ' (l-e-2Asm{ej)sin{ej+,)) (2Asin(^^)sin(0,^+i)) 

Again we can compare this result with our results for generic Nt, expressed in 



section |3.3.1| in terms of character expansions. In the A^^^ = 1 case the sum over 
the representations can be performed exactly, and a closed expression for the Ca/3 
coefficients can be obtained. This can be done by using the identity: 

/(/?)„_! - J(/3)„+i = 2nl{f3)n (80) 

and eq.(|T5D. The result is: 



4 Ji (A) sin(^s) sin(^£+i) 8^/1 (A) sin^ {6^) sin^ y^us+^ j 



8A/i(A) sin2(^,^) sin2(^,.+,) 4Ji(A) sin(^,^) sin(e^ 
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which perfectly matches eg. (|79[) , if we take into account the normahzation by 
and the ^^^^ factor, as in eg. (|76D 

4 Mean Field computation of the critical cou- 
pling 

The effective action obtained in the previous section describes a d dimensional spin 
model with complicated interactions and cannot be solved exactly. However several 
features of the model can be figured out rather easily. First, it can be seen that all 
the interaction terms are even functions of the variables 6', so that the model has 
a global Z2 symmetry, in agreement with the Svetitsky-Yaffe conjecture. Here and 
in the following we shall assume to have fixed the asymmetry ratio p, and we shall 
study the phase diagram of the model in terms of only one coupling (for instance: 
(3t). However, as already anticipated, the p dependence of the critical temperature 
will play a major role. For large values of the coupling j3t the Z2 symmetry is 
spontaneously broken, hence we expect a phase transition for some critical value of 
the coupling. The simplest way to estimate this critical coupling is certainly the 
mean field approximation. As it is well known this method gives in general rather 
rough estimates of the critical temperature, and much more refined technigues have 
been elaborated in these last years. In the section 4.4 we shall comment on this 
point in more detail and apply an improved version of the mean field approximation 
(Bethe approximation) to the A^^^ = 1 case. This test will make us confident of the 
the fact that the errors which we make by keeping a plain mean field approximation 
are of the order of the 10%-15%. Nevertheless, in this section we shall restrict 
ourselves to the plain mean field approximation. The reason is twofold: first, a 
relevant part of our results (and in particular the agreement with the weak coupling 
calculation of Karsch) rely on differences of critical couplings, and these differences 
are only slightly affected by the mean field approximation, which essentially affects 
the data in the sense of giving an overall systematic error. Second, as Nt increases 
the error that we make by using the mean field approximation becomes smaller than 
that due to the truncation of the jSs expansion. The huge amount of complexity 
needed to implement more refined approximations would be justified, and would 
become meaningful, only if higher orders in the expansion were taken into account. 
This could well be done in principle but, as we shall see, the results we obtain by 
keeping only the mean field result are already very interesting, since they clearly 
show the expected trend. Let us also mention that there is another situation in 
which more precise methods to estimate the critical coupling are justified, namely 
the Nt = 1 case, in which the diagrammatic entropy of higher order contributions is 
highly constrained and we can expect that the order alone gives already a very 
good approximation. This will be the subject of sect. 4.4 . 
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4.1 Theoretical expectations. 

The ultimate test of the correctness of any lattice regularization is that, as the 
continuum limit is approached, the various dimensional quantities in which one is 
interested follow the correct scaling behaviour. This scaling behaviour can be easily 
obtained by writing explicitly the dependence on the lattice spacing a of the relevant 
(dimensional) observables. Let us study first the symmetric case (3^ = Pt = P = ^. 
Then the dependence of the lattice spacing on j3 is known in the continuum limit 
in form of the renormalization group equation: 

aAz.= (W)"^exp(^-^^ , (82) 

where is the lattice scale parameter (in units of which we must measure any 
dimensional quantity on the lattice) and bo, bi are the first two coefficients of the 
Callan-Symanzik equation: 

with = 2 in our case. 

Here and in the following we have fixed the spacetime dimensions to be (3+1). 
This is a particularly important remark since it is only in (3+1) dimensions that 
the coupling constant (3 is adimensional and the renormalization group equations 
have this peculiar exponential behaviour. 

Plugging eq.(^) into the definition of critical temperature: 

aNt 

we find: 

Al = iV? (— j ^''P (— j • <«^' 

If the continuum limit is correctly reached then the ratio Tc/A^ should approach 
for large enough values of (3 (hence, in our case, also for large values of Nt) a constant 
value. Plugging this constant into eq. (^5]) we immediately recover the well known 
(approximate) logarithmic growth of (3c as a function of Nt, which is the typical 
signature of the correct continuum limit behaviour of the deconfinement tempera- 
ture in a (3+1) dimensional LGT. It is important at this point to stress that this 
logarithmic behaviour is a non trivial requirement for any effective theory approach 
to the deconfinement transition. For instance, in the zeroth order approximation 
of the effective action, the effective coupling constant is written as a combination 
of modified Bessel functions raised to the Nt power (see eq.(p^). Since the large jd 
asymptotic behaviour of the Bessel function /„(/?) is 
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then the effective couphngs scales as: 



Nt 



[n 



l)Nt 



2(5 



+ 



^7) 



which imphes a linear scahng of f3c as a function of Nt. 

The lack of logarithmic scaling in the zeroth order approximation is one of the 
main reasons which motivated us to look at higher order corrections in the effective 
action. 



4.2 Asymmetric Lattices 

The fact that we have in general asymmetric couplings f3s 7^ A adds some further 
complication to the previous discussion, but has some very important consequences. 
For each value of p we have a new independent regularization scheme, with an 
independent renormalization group equation. This means that, if we define the 
coupling g according to eq.(^ as p- = ^/ PsPt, we must substitute eq.(^2]) with: 

aA(p) = (&o/)"^exp(^-^^ . (88) 

In other words we have that now the A parameter is also a function of p. The last 
step, in order to obtain meaningful results in the continuum limit is then to relate 
the scale A(p) to the A parameter of the symmetric lattice regularization. This 
problem was studied in |^ where the ratio A(p)/Al (with A^ = A(l) denoting the 
scale of the symmetric case Ps = Pt discussed above) was evaluated explicitly at 
one loop. It was found to be a universal function of p: 

where c^- and Cr are the functions already introduced in eq.s (|^ and (|). It is easy 
to see that this result is completely equivalent to the shift in /3 of eq.s (0) and (||). 
It gives us a tool to better understand eq.(|^), which simply encodes the effect of 
the quantum fluctuations at one loop. 

This result is particularly interesting for our purposes, since it allows to extend 
our analysis also to asymmetric lattices, taking into account also the quantum 
effects. Let us discuss in detail this point. 

Consider an asymmetric regularization, with /3t = p'^Ps, P > 1, defined on a 
lattice with temporal extension Nt. By using eq.(|D and (H) we see that this regu- 
larization is (classically) equivalent to that on a symmetric lattice with (3 = ^ = (3gp 

and with temporal extension Nt = Hence p must be a rational number; we shall 
choose in general p to be an integer number, so that A^^ will be a multiple of Nt. 
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At the quantum level we can still show the equivalence of the symmetric and asym- 
metric lattice regularizations provided we modify the previous relations according 
to eq.(|^ and (||). By using the explicit knowledge of the p — > cxd expansion of 
and Cri we obtain, for large enough values of p, the following scaling behaviour: 

f3t,c{p) = {(^c + oP,)p + a\ (90) 

where Pt,c{p) is the critical coupling on the asymmetric lattice and (3c the critical 
coupling on the equivalent symmetric lattice. The numerical values of and al 
are reported in sect. 2, following eq.(^. Higher order corrections to eq.(|90D vanish 
as p — > CX3. 

The limit p ^ cxd is particularly interesting because some relevant simplifications 
occur in that limit in the effective action. Let us consider first the contributions 
of the timelike plaquettes in the zeroth order approximation given in eq.(^^. We 
have now to perform in (|2^) the following substitution: 



Nt^Nt = pNt , Pt^Ptip) (91) 
with (3t{p) given in eq. (|90|) , and finally take the limit p — oo . By using the 



asymptotic expansion of the Bessel functions (|86D one obtains in that limit that the 



effective couplings at the r.h.s. of eq.(p^) become 



him) 



p^oo , 2j(j + l)iV, 
- exp(— ^^-^), (92) 



where f3 is the coupling of the corresponding symmetric lattice. One can easily 
recognize the quadratic Casimir of the representation j at the exponent, and one 
immediately realizes that the action for the time-like plaquettes becomes the Heat- 
Kernel action in the p ^ oo limit, namely in the hamiltonian limit. 

The couplings of the spacelike plaquette in the adjoint representation and of the 
pair of plaquettes in the fundamental representation are respectively proportional 



to SNt ^/fA^l] and AN^ C/ff]] )^ (see eq.(^)). We replace A^^ as in (^) and /?,(p) 



by and take the limit p — > oo. The contribution of the adjoint representation 
is of order 1/p and hence it vanishes in this limit, due essentially to its zero measure 
in an infinite lattice. On the contrary the coupling of the pair of plaquettes in the 
fundamental representation is finite in the limit p — oo and its limiting value is 
Nf{P + a°)^/4. Notice that all powers of Psip) higher than Ps{pY in the power 
expansion of the Bessel functions give a vanishing contribution in the limit p ^ oo. 
This means that, as already remarked at the beginning of sect. 3, the character 
expansion and the power expansion coincide in this limit, thus removing any possible 
ambiguity. Notice also that in spite of the vanishing of /?s(p) in the asymmetric limit 
the effective expansion parameter Nt{P + a°) is never very small and it does indeed 
increase as we approach the continuum limit. 
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p 


Nt 






6 


7 


1 


6 


3.216 


3.216 






2 


12 


5.882 


1.471 


2.666 


0.55 


3 


18 


8.576 


0.953 


2.694 


0.49 


4 


24 


11.282 


0.705 


2.706 


0.46 


5 


30 


13.993 


0.560 


2.711 


0.45 


6 


36 


16.707 


0.464 


2.714 


0.42 


7 


42 


19.422 


0.396 


2.715 


0.42 


8 


48 


22.138 


0.346 


2.716 


0.41 


10 


60 


27.571 


0.276 


2.717 


0.39 


20 


120 


54.753 


0.137 


2.718 


0.39 



Tab. I The critical coupling (3t,c as a function of p. In the second column we 
have reported the temporal extension Nt = pNt of the asymmetric lattices that we 
used. In the fourth column we have reported the values of f3s^c = A,c/p^- In the 
last two columns we have reported the values of'-/ and 6 (see text for explanation). 
The values of 6 reported in Tab. I are obtained by applying the definition / jg^ j with 
p2 = p o,nd pi equal to the value of p in the previous row. 

4.3 Results. 

In order to extract reliable estimates for the critical couplings, we must first discuss 
the behaviour of our results as functions of the asymmetry parameter p. 

4.3.1 p dependence of the results 

It is very interesting to study the p dependence of our results, to see if eq.(|90D 
is fulfilled. In Tab. I we have reported, as an example, the p behaviour of a set 
of asymmetric regularizations which are all equivalent to the symmetric Nt = 6 
case. For each pair of (subsequent) values of Pt,c{p) we have constructed the two 
quantities 

P2- pi 

7 = A,c(p2) - Sp2 (94) 

which, if eq. (|90|) is fulfilled, should provide a good estimate of Pc + and a] 
respectively. Their values together with those of (3t,c are reported in Tab. I. 

It can be seen that the data follow very well the expected law. In particular it 
is clearly visible the 1/p quantum correction, which is definitely different from zero 
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Nt 


5 


e 


7 


2 


1.554 


-0.184 


0.414 


3 


1.971 


-0.210 


0.375 


4 


2.259 


-0.221 


0.373 


5 


2.500 


-0.235 


0.372 


6 


2.718 


-0.249 


0.389 


8 


3.114 


-0.271 


0.413 


16 


4.443 


-0.327 


0.508 






-0.27192 


0.50 



Tab. II Values of 6, e and j as functions of Nt (see the text for the definitions 
of these three quantities), e and 7 are estimators of and a\ whose theoretical 
values are reported, for comparison, in the last row of the table. 

and smoothly approaches for large values of p the value 7 ~ 0.39, which is not too 
far from the expected value a\ = 0.5. 

The values of j3s,c = Pt,c/p'^, that we have reported in the fourth column of Tab. 
I, give an idea of the reliability of our strong coupling expansion in /J^. 

We repeated the same analysis for all the values of Nf for which Montecarlo 
data arc known (the MC data are reported in Tab. IV). Our results are collected 
in Tab. 11 where we have reported the asymptotic (large p) values of 7 and S as 
well as of another quantity e that should provide us with an estimate of a^. This is 
defined as e — 5 — (3c, where (3c is the critical coupling obtained with a symmetric 
lattice of size Nt by keeping strictly only the terms of order in the character 
expansion. In other words, in computing f3s we neglect all contributions that vanish 
in the limit p 00, thus keeping the same contributions in the symmetric and 
asymmetric lattice. Tab. II shows that, in the range Nt = 2 — 16, the agreement 
with the theoretically expected values of q;° and a] is really remarkable. Let us 
stress again that this agreement is highly non trivial since q;° and were obtained 
with a weak coupling calculation, while our effective action is the result of a strong 
coupling expansion. The reason of this success is very likely related to the fact that 
we have been able to sum to all orders in (3t the timelike contribution of the effective 
action. 

4.3.2 Scaling behaviour. 

The agreement between the p dependence of our mean field results and the the- 
oretical expectations allows us to be confident on their consistency, also at the 
quantum level. So in order to extract our best estimate for the critical coupling we 
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Nt 

Fig. 2 Values of the critical coupling f3c are plotted for different values of the 
number of time-like links Nt. Results obtained with Montecarlo simulations, which 
are denoted by *, are compared with those obtained with our mean field analysis: A 
represents the data for Pdo, the critical coupling in the zeroth order approximation 
and □ the data for f3c\i,p^oo, the critical coupling including the effect of the space- 
like plaquettes al the lowest non trivial order, calculated in the limit p — > oo and 
reduced to the value p = 1 as described in the text. 

can consistently assume eq.s (|^) and witli tlie correct theoretical a coefficients, 
work in the large p limit (so as to avoid ambiguities in the expansion) and then 
rescale the resulting critical values to the limit of symmetric lattice^, so as to allow 
a comparison with the Montecarlo results. 

The values obtained in this way for the critical couplings are reported in Tab. 
Ill, and plotted in Fig. 2, where they are also compared with the Montecarlo results 
(extracted from which are reported in Tab. IV . It is impressive to see how the 
contribution of the space-like plaquettes, although taken into account at the order 

only, improves the agreement of the results of the mean field method with the 
ones of the Montecarlo simulations. To the order (3^ the mean field method gives 
results which are in reasonably good agreement with the Montecarlo for A^^t's as high 
as 5, displaying in the range 1 < iVj < 5 a behaviour which is compatible with the 
logarithmic scaling predicted by the renormalization group. In contrast the model 
with only time-like plaquettes shows from A'^t = 1 a linear scaling behaviour in con- 

^Let us notice, as a side remark, that the difference between the values of /?c obtained in 
this way and those which one would obtain with the naive procedure of choosing right from the 
beginning a symmetric lattice is not very large, but nevertheless it is not negligible and could 
become very important if further orders in the Ps expansion were to be added to the effective 
action. 
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Nt 


Pc\o 


Pc 1 l,p— »oo 


2 


1.957 


1.723 


3 


2.957 


2.089 


4 


3.876 


2.366 


5 


4.763 


2.606 


6 


5.636 


2.826 


8 


7.369 


3.227 


16 


14.262 


4.572 



Tab. Ill The critical coupling jSc as a function of the lattice size Nt in the t 
direction. In the second column the values obtained with the zeroth order approxi- 
mation and in the third column those obtained at the order (3l by rescaling the values 
obtained in highly asymmetric lattices p ^ oo as explained in the text. 

tradiction with both renormahzation group predictions and Montecarlo simulations. 
It must be stressed however that the almost perfect agreement of our results with 
the Montecarlo simulations for small A^^^'s (especially for Nt in the range (3 — 4) ) is 
most probably only apparent, resulting from a compensation between an expected 
10-15% systematic error due to the mean field approximation and the effects of 
higher orders contributions. 

4.4 Nf = 1 again: improving the mean field approximation. 

As we remarked in sect. 3.5, the most interesting feature of the Nt = 1 case is 
that, due to the remarkable simplifications which occur in this case, the character 
expansion of the effective action can be resummed exactly. This allows a much 
simpler implementation of improved versions of the mean field approximation, and 
we decided to use this laboratory to test these improved estimators and 

above all to have a hint of the magnitude of the systematic errors involved in the 
plain mean field approximation that we used in the previous section. Notice also 
that for this particular Nt = 1 case it does exist a very precise Montecarlo estimate 
of the critical coupling, namely Pc = 0.8730(2) |l^, and we shall use this value as 
a reference point to compare our predictions. 

Let us notice, as a preliminary remark, that in the Nt = 1 case the contri- 
bution is given by the adjoint plaquette term only, since there is no room to locate 
a pair of plaquettes in the fundamental representation. This fact has two conse- 
quences: first, the magnitude of the correction to the critical temperature due to 
the spacelike contribution is much smaller than in the A^^^ > 1 case; second, it has 
the opposite sign, namely (3t,c increases as a consequence of the spacelike term. 
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S 

He 


T /A 


2 


1 88nn('3n') 


29 7(2) 


3 


2.1768(30) 


41.4(3) 


4 


2.2986(6) 


42.1(1) 


5 


2.3726(45) 


40.6(5) 


6 


2.4265(30) 


38.7(3) 


8 


2.5115(40) 


36.0(4) 


16 


2.7395(100) 


32.0(8) 



Tab. IV The critical coupling f3c and the corresponding deconfinement temperature 
Tc/A as a function of the lattice size in the t direction, Nf, in the (3+1) dimensional 



SU(2) LGT. The data are taken from / (7^/ . 



The simplest way to improve the mean field approximation is to consider larger 
and larger clusters of spins (see Fig. 3). It can be shown that this modification 
indeed allows a more and more precise determination of the critical coupling. The 
first improvement (step 2 in the notation of Fig. 3) is also known as Bethe approx- 
imation |T^. In Tab. V we report the results of our analysis. It is easy to see that 



both with and without the spacelike plaquettes the plain mean field approximation 
is affected by an error which ranges from 10% to 15% (depending on the role which 
is played by the spacelike contribution in the Montecarlo estimate), and that almost 
half of this gap is filled by the Bethe approximation. 

We will show in a forthcoming publication ||18[, that the results of next order 



approximation (step 3 in the notation of Fig. 3) differ from the exact result by only 
a few per cent. Notice however that the step 3 approximation is very hard to handle 
and requires several technical manipulations. In this forthcoming paper, we will also 
show that a result which is very close to the Montecarlo one given in last column of 
Tab. V can be obtained by comparing the different steps of Bethe approximations 
in the present model with the corresponding ones in the Ising model, whose critical 
coupling is known with high precision. 

An alternative approach one can follow is to identify explicitly the Ising model 
which is hidden in the effective action (following the Svetitsky-Yaffe analysis) and 
then use again our knowledge of the critical coupling of the three dimensional Ising 
model. This approach has been developed for the Nt = 1 case in M and it will also 



be fully exploited in [18 
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Fig. 3 Progressively more precise Bethe-like approximations 



step 1 step 2 MonteCarlo 

0.7702 0.8139 

0.7705 0.8145 0.8730(2) 



Tab. V The critical coupling (3c for the action Sq (first row) and Sq + 5*1 (second 
row) in the mean field approximation (step 1) and the Bethe approximation (step 
2) . In the last column the Montecarlo result of ref. (including full contribution 
of space-like plaquettes). 

5 Conclusions 

The main result of our paper is the construction of the effective action for the 
Polyakov loops at the first non trivial order in the spacelike coupling, which is 
exact to all orders in the timelike coupling. 

Extracting the critical temperature from this effective action, and comparing 
our results with those of the Montecarlo simulation, we have seen that the effective 
action describes rather well (even within the mean field approximation) the full 
theory up to a lattice size in the compactified time direction of iV^ ~ 5. This is an 
impressive improvement with respect to the previous studies, in which the spacelike 
contributions were neglected and which were constrained to Nt = 1. We see in 
principle no obstruction to extend our analysis to higher orders and to reach a 
better and better agreement with the renormalization group expectations. Such a 
result would obviously be very interesting, since it would be the first time that one 
can reach the scaling regime of a dimensional quantity (as the critical temperature) 
by using only analytical tools. From this point of view the critical temperature 
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seems indeed in a better position with respect for instance to the string tension, for 
which, due to the roughening transition, a strong couphng approach has very few 
chances to reach the scahng region. It is also remarkable the agreement between 
our results, based essentially on a strong coupling expansion, and the weak coupling 
calculations of Karsch, regarding the consistency of the theory with respect to lattice 
deformations described by the asymmetry parameter p. 

Moreover our approach can be straightforwardly extended to the SU(3) gauge 
model, which is obviously more interesting from a phenomenological point of view. 

Besides these obvious remarks there are two more reasons of interest of our 



result, which we think should deserve some attention. First, it was noticed in 
that the ratio T^/ ^/a shows a very precocious scaling, and is essentially stable 
already for lattices as small as Nt = 4. This value lies inside the region where we 
have seen we can trust our expansion. Therefore, if one were able to extract the 
string tension in our framework, it would be possible to reach the continuum limit 
value of Tc/ ^/a already within our (relatively simple) effective action. 

A second possible interesting application is in the study of LGT with suitable 
generalizations of the Wilson action. In particular, much interest has been recently 
devoted to the so called fundamental-adjoint SU(2) LGT (see [0 and references 
therein) in which a term proportional to the plaquette in the adjoint representation 
is added to the ordinary Wilson action. This introduces a new coupling Pa- In 



the behaviour of the deconfinement transition in the extended coupling plane (/?, (3a) 
was studied with a strong coupling effective action truncated to the first order and 
for Nt = 2 and Nt = 4, hence in a region in which our approach seems to have a 
good behaviour. Since our action is written in terms of a character expansion, it 
should be rather straightforward to generalize it to include a coupling Pa different 
from zero. 
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